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*-^H, I give an account of a new definition of transverse-momentum-dependent parton densities. 

*-». , The new definition solves a number of difficulties and inconsistencies in earlier definitions. 

<D ' 

1. Introduction 

In this talk, I presented a new definition of transverse-momentum-dependent (TMD) 
parton densities that I developed in my recent bookM-' They overcome many defi- 
ciencies of previous definitions. Numerical values for these densities with the new 
definition were obtained by Aybat and Rogers^ from fits to data. 

Before embarking on the QCD definition, I give a reminder about how TMD 
[ — ■ parton densities arise naturally in the parton model. Then I summarize some of 

the difficulties that have to be overcome in defining TMD parton densities in QCD. 
After that I present the new definition, compare them with other work, and indicate 
some of the implications of the new definition. 

A number of assertions are given here without justification. For more details on 
the justification, see Ref. Q] 



2. Basics of TMD parton densities 

A motivation that TMD parton densities are useful objects to use is provided by the 
parton model for the Drell-Yan process, i.e., the production of a high-mass lepton 
pair in a high-energy hadron-hadron collision. In the parton model, the lepton pair 
is formed by annihilation of a quark in one hadron and an antiquark from the other, 
Fig. [Q This gives TMD factorization: 

■jjAj^ = Yl J d2k AT f j/hA (x A ,k AT ) fj /hB (x B ,q T -k AT ) -^-, (1) 
3 

where the fj/H factors are the TMD parton densities, functions of a longitudinal 
momentum fraction and a transverse momentum, while da is the partonic hard 
scattering in the lowest-order approximation. The query over the equality sign is a 
reminder that this is not exactly the correct TMD factorization theorem in QCD. 
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Fig. 1. The parton model for the Drell-Yan process 

Because the transverse momentum q T of the lepton pair is the sum of the trans- 
verse momenta of the partons, the cross section (JlJ is directly sensitive to partonic 
transverse momenta. If the partons had no transverse momentum, the cross section 
would be a delta-function at q T = 0. This contrasts with deep-inelastic scattering 
where only one parton participates in the hard scattering, so that its transverse 
momentum can be neglected with respect to the large momentum transfer Q in the 
hard scattering. 

The need for TMD parton densities in Eq. ([1]) establishes that TMD parton 
densities are important quantities for a quantitative description of many hard pro- 
cesses. 

To derive the parton model, one needs to use a cancellation of spectator- 
spectator interactions!^ In addition one needs to assume other topologies of graph 
are unimportant, that partonic fcp and virtuality are limited, and that no higher- 
order corrections are needed to the hard scattering. All of the last three assumptions 
are violated in QCD and are associated with a need to modify the definitions of the 
parton densities and the factorization formula in QCD. 



2.1. Explicit definition of TMD parton density: complications in 
QCD 

In constructing an operator definition of a TMD parton density in a hadron, I 
assume that the hadron is moving in the +z direction, and I will use light-front 



coordinates defined by v^ — (v^ 



, v T ), with v ± = (v° ± V z )/y/2, v T = (v x ,v y ). 



The parton model leads to a definition of a parton density as a hadron expecta- 
tion value of the number density of a parton, as specified in light-front quantization. 
A first attempt at applying this in QCD uses the A + = gauge. This is equiva- 
lent to the following gauge-invariant definition with a Wilson line in the direction^ 
-n = -(0,l,0r): 



a Note that the derivation of factorization requires that parton densities for the Drell-Yan process 
use past-pointing Wilson linesptl 
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Fig. 2. Wilson line in naive definition of gauge invariant parton density, viewed from the side. 



fj/h(£,kr) 



(2^) 4 * 2 




= F.T. (P\ i> 3 (0,w',w T ) W[w, -n] 
where the Wilson line is 



t r 



[tr.link]W[0,-7#,-(0)|P) , (2) 



W(co, x; —n) = P exp 



-*5o 



ds (-n) • Ag(x-sn)i a 



(3) 



"F.T." denotes a Fourier transform with respect tom" and i«t, and "[tr. link]" 
denotes a transverse Wilson line at infinity. The product of Wilson lines is a path- 
ordered exponential of gluon fields along the line shown in Fig. [2] 

There are several complications that prevent this definition from being taken 
literally: 

• UV divergences. These are removed by suitable renormalization countert- 
erms. 

• Rapidity divergences!^ These concern gluons whose rapidity goes to — oo. 

• There are Wilson line self energies. These are present in the above definition 
of the TMD parton density, but not in the actual cross section. They create 
divergences as the length of the Wilson line goes to infinity. 



2.2. Example: Regions for one gluon 

These complications are illustrated by the graph in Fig. [3f a), where into a parton- 
model graph for the Drell-Yan process is inserted one gluon exchanged between 
the annihilating quark and antiquark. Associated with this in the factorization of 
collinear configurations is Fig. [3jb) for the parton density in hadron A and a similar 
graph for the parton density in hadron B. 

These graphs give leading power contributions from regions of gluon momentum 
indicated by the shaded/blue regions in Fig. 01 a plot of the space of gluon rapidity 
and the logarithm of its transverse momentum. A frame is used in which hadron A 
moves in the +z direction and hadron B in the — z direction. The red line indicates 
roughly the region needed when the transverse momentum q^ of the lepton pair is 
fixed. The leading logarithm approximation is obtained by replacing the integrand 
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(a) 

Fig. 3. Graphs for the Drell-Yan cross section and for a TMD parton density in hadron A. 






Fig. 4. Regions of gluon momentum dominating in the graphs of Fig.[3]and in the corresponding 
graph for the TMD parton density in hadron B. 



in J d In h? dy . . . by a constant. In the cross section integrated over qx, the trian- 
gular region in the first plot gives the well known Sudakov double logarithm of Q 
for Fig. \%l&) . The edges of the triangle are set roughly by the rapidities of the an- 
nihilating quark and antiquark, and by an infra-red cutoff on gluon fcx- In massless 
perturbation theory, with m g = 0, an extra doubly logarithmic contribution comes 
from even smaller fcx- 

Replacing the gluon and the lower part of the graph by a corresponding graph 
for the parton density in hadron A, Fig. (21b), gives the region shown in Fig. EJb). 
The parton density graph corresponds to a good approximation for the gluon in 
Fig. [3ja) when it is collinear to hadron A, i.e., when y is sufficiently positive and 
fcx is sufficiently much less than Q. The proof of this assertion uses a Ward identity 
to move the upper end of the gluon line in Fig. EJa) onto a Wilson line. 

But the approximation becomes very bad for negative y, so much so, that there 
is a divergence in the integral to y = — oo (a "rapidity divergence"). 

Similarly, in Fig. SJc), the corresponding gluonic term in the parton density in 
hadron B reproduces the region of gluon momentum collinear to hadron B, but 
fails badly in the positive rapidity region. 

In the QCD factorization theorem that generalizes Eq. (JT]), these two parton- 
density contributions are to be added. A further problem is then that the region of 
central gluon rapidity is double counted. 

All of the problems are overcome by finding a suitably modified definition of the 
TMD parton densities, together with a valid factorization theorem using the parton 
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densities. 



3. TMD factorization in QCD 

TMD factorization in QCD has the general form given by Collins, Soper, and 
StermarPI (CSS): 

da = Hx convolution of ABS + high-gx correction (Y) + power-suppressed. (4) 

Here, if is a hard-scattering factor, A and B are TMD parton densities, and S is 
a soft factor, with the product ABS being a convolution in transverse momentum. 
Double counting subtractions are used in the definition of the factors, and there 
must be present some kind of cut off on rapidity divergences. 

The HABS part correctly gives the cross section when qx <§C Q, but it fails 
at larger <yr- The Y term corrects the errors at large qx; it has the form given by 
ordinary collinear factorization, but with a subtraction of the low gx region. 

All errors in Eq. Q are suppressed by a power of 1/Q. The pattern of errors in 
the individual terms is interesting: 



Errors in 



collinear factorization 

TMD factorization (HABS) 

TMD factorization + Y 



are a power of < 



A 

q T 

A q T 



Q' 

Each of the A, B, and S factors has a non-perturbative contribution that cannot 
currently be predicted from the theory; to determine the non-perturbative contri- 
butions one must fit them to data. Unfortunately, the soft factor S always appears 
multiplied by two collinear factors (not only in the Drell-Yan process, but also in 
similar factorization theorems for other reactions). Thus there is no possibility of 
measuring S independently of the parton density factors A and B. 

A better formulation can therefore be obtained by redefining the parton density 
factors so that they each incorporate a square root of the soft factor. The new 
definition, presented below, accomplish this, albeit somewhat indirectly. 



4. New definition 

The new definition is in terms of unsubtracted quantities that have non-light-like 
Wilson lines. The unsubtracted parton density in hadron A has a Wilson line of 
rapidity y 2 : 



runsub 
Jf/H A 



(x,b T ;yp A -y 2 ) = Tr Tr — 

color Dirac Z 



+ f dk~ d 2 ~ 2e k 



(27T) 



■1 „ — iktT-bT 




6 Collins 



The unsubtracted soft factor has Wilson lines of rapidities y\ and 2/2 j with y\ > y 2 ; 



SOr)*±- r ' r 



N C J (2tt) 



KS e -ik S T-bi 



4-2e 




(7) 



The definition of the TMD parton density uses an auxiliary rapidity parameter y n 
which sets the finite rapidity of certain Wilson lines. The rapidities of other Wilson 
lines are taken to +00 or —00: 



clef 



ff/H A (x, b T ; Ca; aO =' „,. lim_ ff/gZ(x, b T ; VPa -y 2 ) 



S(b T ,yi,y n ) 



3/1-H-oo 
8/2- 



S(b T ,y 1 ,y 2 ) S (0) (b T ,y n ,y 2 ) 



S(b T ; +00, y n ) 



/— (x,6 t;2/pa -(-oo)),U— 

y ,S(&t;+oo, -00) A( )(0T;yn,-oo) 



(8) 



Here ( A is defined to be M 2 PA x 2 e 2< - VP A-v^ to match a definition in Ref. 1 A UV 
renormalization factor is implicit, as is the limit e — >• for the removal of the 
dimensional regularization of UV divergences. 

The multiplications, divisions, and square roots are applied to the various factors 
in transverse coordinate space. They would correspond to complicated convolutions 
were they to be applied in transverse momentum space. 

Diagrammatically, this definition can be written as 



(9) 





4.1. Why this strange definition? 

This definition seems unexpectedly complicated, although the square root is to be 
expected given the motivation of absorbing a square root of the soft factor into 
each parton density in the factorization property. In fact, the definition is unique 
(up to the choice of the prescription for renormalization of UV divergences) given 
the following properties. 

• The parton density (in transverse coordinate space) is the product of the 
basic parton density and powers of the soft factor. 

• Light-like Wilson lines are used if possible: 

— The basic parton density has only a light-like Wilson line 
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— Each soft factor has at most one non-light-like line 

• Rapidity divergences cancel, as do divergences when the length of the Wil- 
son lines goes to infinity (L — > oo). 

Note that only factorization-compatible definitions should be considered, that the 
above definition is factorization-compatible, and that the definition implements 
double-counting subtractions between the soft region and the two collinear regions 
of momentum space. 

4.2. Consequences of the strange definition 

Some consequences of the new definition are: 

• Factorization holds, as in Eq. (|4]), but without soft factor. 

• Thus the TMD part of factorization is just like that in the parton model, 
{l}, except that the parton densities depend on C,a (and a corresponding 
variable £b for the parton density in hadron B) and on (i, and except that 
the hard scattering has higher order corrections. Note that OiCb — Q 4 - 

• Rapidity divergences and Wilson-line self-energy divergences cancel. 

• Contributions from the gauge link at infinity cancel in Feynman gauge. 

• There is an effective cutoff on gluon rapidity at y n 

• When the energy of an experiment and the dilepton mass are scaled up, 
the parton densities change because (a and Cs increase. The dependence 
on these variables is governed by the dependence on y n , the rapidity of the 
non-light-like Wilson lines, and this dependence is governed by the CSS 
evolution equations to be discussed below. 

In addition, the renormalization scale u should be of order Q to avoid 
large logarithmic higher-order corrections in the hard scattering. The fj,- 
dependence is governed by a renormalization-group equation. 

• The versions of the CSS evolution equations for the newly defined parton 

[7l 

densities are simpler than the original ones. In particular, they no longer 

have hard-to-control power-law corrections. 

5. Evolution, etc 

In this section, I summarize the evolution equations and the small-frx expansions. 
These enable the TMD parton densities to be expressed in terms of (a) non- 
perturbative quantities without scale dependence, and (b) perturbative quantities 
that do not have large logarithms in the coefficients of their perturbative expansions. 
The dependence of the parton density $8§ on (a is obtained by differentiating 
with respect to y n the two soft factors that use a Wilson line of rapidity y n . This 
gives a CSS-style equation 

^^y ^ ^^ (io) 



8 Collins 

The renormalization group (RG) equation of the kernel of this equation is 

dK , , . . , 

^ = -7*(. 9 (.)). (ID 

The RG equation of the TMD parton density is 

-^ =7(,s(w;C/m )• (12) 

From these equations, follows an equation for the ^-dependence of 7, and hence that 



t(<?(m); C/m 2 ) = 7(5(m); l) - ^(<?(m)) in -V- (13) 



At small 6t, the TMD parton density has a factorized formula in terms of ordinary 
integrated parton densities: 

///#(«, &t; C; m) = 5Z / —Cf/j(v/x,bT;C,fi,g(p)) fj/H&v) + o[(m6 T ) p ] 

See Aybat's talk at this conference for results at level of TMD parton densities. 



6. Comparisons 

In this section, I compare the new definition with some other definitions that have 
appeared in the literature. 

6.1. CSS- style definitions 

The new definition can be considered as a development of that originally proposed 
by Collins and SoperPThe old definition used the basic operator definition @, but 
without a Wilson line and with the use of a non-light-like axial gauge n ■ A = to 
cut off rapidity divergences. 

With this definition there was a lack of an actual proof of TMD factorization 
for the Drell-Yan processM^ The problem is that Glauber region is difficult to treat 
in "physical gauges" such as the n ■ A = gauge. This lack has been remedied in 
Ref. [I] with the aid of the new definition. The old definition can be converted to 
a gauge-invariant form, with the use of non-light-like Wilson lines, which should 
be space-like to get factorization!^ However there are divergences associated with 
self energies on the infinitely long dipolar Wilson lines. This last problem was only 
recently noticed, by Bacchetta, Boer, Diehl, and Mulders — see App. A of Ref. 6. 

Furthermore, with the Collins-Soper definition there was a separate soft factor 
in the factorization formula, but as stated earlier, the non-perturbative part of the 
soft factor cannot be independently determined from data. 

Finally the derivation of the version of the evolution equation ()10|) appropriate 
to the old definition involved ignoring errors that are power-suppressed only at small 
transverse momentum. Since the TMD densities are also used at large transverse 
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momentum (in conjunction with the Y term in (0|), it follows that the actually- 
used TMD pdfs in the Collins-Soper method are not actually those that Collins and 
Soper defined explicitly. 

Ji, Ma, and YuarP improved the Collins-Soper definition of the TMD densities 
and the soft factor with the use of multiple non-light-like Wilson lines. However, 
they did not see how to take infinite-rapidity limits as in Eq. (j5]), which would have 
simplified their formalism. Their factorization formula still contains a soft factor. 



6.2. Naive light-cone-gauge definition 

A number of authors have defined a TMD density as the naive expectation value of 
a light-front parton number operator in light-cone gauge, i.e., by Eq. ([2]) in A + = 
gauge with omission of the Wilson line. Such a definition suffers from rapidity 
divergences P It can only be considered a valid definition if a cutoff is applied. But, 
as far as I can see, the need for a cutoff is often not recognized, or, if it is, a cutoff 
is applied in such an implicit way that it is difficult to work out what the actual 
definition is. 

Similar problems afflict much work that uses TMD densities at small x with 
BFKL-related approaches — see Avsar's talk for more details. 



6.3. Cherednikov and Stefanis's work 

In several recent papers, Cherednikov and Stefani a ' ' ' have presented an anal- 
ysis of possible definitions of TMD parton densities. They have drawn conclusions 
about the anomalous dimensions of the TMD parton densities and about the role of 
gauge links at infinity. Particular results concern the presence or absence of rapid- 
ity divergences in certain versions of light- cone- gauge quantization of QCD. Their 
results are based primarily on one-loop calcu lations. In particle, they claim that 
the Mandelstam-Leibbrandt prescriptio n I I is particularly suitable, by giving fi- 
nite results in situations where other methods of using the light-cone gauge give 
divergences. 

Their gauge-invariant definition is (slightly modified from) 





fffgf^ixte) 



(14) 
with a certain treatment of gauge links at infinity. In calculations light-cone gauge 
is used, with a particular cutoff in the gluon propagator: 



^ML 2 



ff" + (q^n" + n^q") 



1 



1 



2 \q 



vq 



ii] 



(15) 
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where 77 nonzero and positive. This definition in fact has uncanceled rapidity diver- 
gences from graphs that the authors appear not to calculate. Examples are 



(16) 

/ 

(Note that the regulated propagator (|T5|) does not obey ^flj^ = 0, so in these 
graphs the gluon attachments to the Wilson line of negative infinite rapidity are 
nonzero.) Furthermore Cherednikov and Stefanis do not analyze how their definition 
is to be used in a valid factorization formula. 

6.4. Definitions in the soft-collinear effective theory (SCET) 
framework 

In work in SCET, a number of authors have tried to make definitions of TMD 
functions. There are two approaches. 

One is that of Becher and NeubertM^They use Smirnov's subtraction methods^ 
applied in a style appropriate for SCET. They are able to define the product of TMD 
parton densities used in the QCD TMD factorization formula, i.e., they define the 
product ABS in Eq. (J4J . Although the product is called the product of two TMD 
parton densities, there are, in fact, subtractions that effectively give a soft factor. 
But Becher and Neubert are unable to define individual TMD parton densities. As 
far as I can tell, their product of TMD parton densities is the same as with my new 
definition; I think this is guaranteed by the use of Smirnov's subtraction methods. 

The other school in the SCET community is represented by Mantry and 
Petriello! II They work with what they call an "impact parameter beam func- 
tion" (iBF). Compared with the definitions such as ([2]) and ((6]), the fields in the 
definition of an iBF are not restricted to a null plane. So the iBFs are actually more 
general objects than TMD parton densities, and in some other literature would be 
called fully unintegrated parton densities. There are processes more exclusive than 
the Drell-Yan process for which their use is advantageous or even necessary. The 
iBFs have zero bin subtractions, which are something like the soft factors in my 
definition (jSJ. However, the exact relation is not clear to me. 

7. Outlook: Implications of the new definition and associated 
derivations 

• We now have a precise satisfactory operator definition of TMD parton den- 
sities that can be taken literally. 

• The gauge links at infinity are no longer needed in Feynman gauge. 

• The definition is part of a new formalism with much better subtraction 
methods^ than were previously available. 
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There is now a full proofi- of TMD factorization for the Drell-Yan process. 

(CSS^gave up on this!) 

Hence, there is now an fully unambiguous method for calculation of the 

hard scattering coefficient. 

The formalism is much cleaner mathematically than previous versions. 

Although I did not explain this, there is a fully specified relation between 

exact and approximated parton kinematics in the hard scattering. 

As for future work, there is an urgent need to numerically relate TMD 

parton densities with the new definition to those in other formalisms. 
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